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ABSTRACT 

We present phase resolved optical spectroscopy and X-ray timing of the neutron star 
X-ray binary EXO 0748—676 after the source returned to quiescence in the fall of 2008. 
The X-ray light curve displays eclipses consistent in orbital period, orbital phase and 
duration with the predictions and measurements before the return to quiescence. Ha 
and He I emission lines are present in the optical spectra and show the signature of 
the orbit of the binary companion, placing a lower limit on the radial velocity semi- 
amplitude of K% > 405 km s -1 . Both the flux in the continuum and the emission lines 
show orbital modulations, indicating that we observe the hemisphere of the binary 
companion that is being irradiated by the neutron star. Effects due to this irradiation 
preclude a direct measurement of the radial velocity semi-amplitude of the binary 
companion; in fact no stellar absorption lines are seen in the spectrum. Nevertheless, 
our observations place a stringent lower limit on the neutron star mass of Mi > 
1.27M . For the canonical neutron star mass of Mi = 1.4Mq, the mass ratio is 
constrained to 0.075 < q < 0.105. 
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1 INTRODUCTION 

Neutron star low-mass X-ray binaries (LMXBs) consist of a 
neutron star that accretes matter from a low-mass compan- 
ion star. From binary evolution models, we know that the 
neutron star could accrete up to 0.7 Mq from the compan- 
ion (e.g. van den Heuvel & Bitzaraki 19951. In that case, 
neutron stars as massive as > 2 Mq may be expected in 
LMXBs. Such a high-mass neutron star would pose strong 
constraints on the equation of state (EOS) of matter at the 
pressures and densities encountered in neutron star cores 
(see Lattimer & Prakash 2004 and references therein). De- 



termining the neutron star EOS is one of the ultimate goals 
of the study of neutron stars. 

In order to determine the neutron star mass in these ef- 
fectively single-line spectroscopic binaries, we need to know 
the inclination, the projected rotational velocity vsini of 
the companion star and the semi-amplitude K2 of its radial 
velocity curve. In principle, these last two parameters can 
be determined via optical spectroscopy. Unfortunately, the 
companion stars of LMXBs are nearly always intrinsically 
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faint and they are typically located at distances of several 
kilo-parsecs. Furthermore, in systems accreting at high ac- 
cretion rates the accretion disc out-shines the small compan- 
ion stars making it impossible to detect the stellar absorp- 
tion lines necessary to determine K2 and vsini. In LMXBs 
the binary inclination can be determined accurately if the 
system is eclipsing (Hornc 19851. Hence, in order to provide 
the strongest, model independent constraints on the neutron 
star mass we need to target quiescent, eclipsing LMXBs. 

EXO 0748-676 is such an eclipsing LMXB. It was dis- 
covered using observations of the European X-Ray Observa- 
tory Satellite (EXOSAT) in 1985 ( |Parmar et al.|1985[ ). The 
presence of X-ray eclipses and type I X-ray bursts marking 
the compact object as a neutron star, were readily found 
(Parmar et al. 1986 Gottwald et al. 19861. In hindsight 



the source had been observed serendipitously by the EIN- 



STEIN satellite as well as during EXOSA T slews ( |Garcia fc 
Callanan|1999 and Reynolds et al.|["l999 respectively). 



The optical counterpart to the X-ray source has been 



found by Parmar et al. (19851. Wade et al. (19851 report 



on the discovery of the optical counterpart as well as on 
searches of archival optical plates showing that the optical 
source was not present down to roughly 23rd magnitude in 
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the SRC J plate when the X-ray source was inactive. Cot- 
( 2002 I reported on a measurement of the grav- 



tam et al. 



itational redshift from X-ray absorption line spectroscopy 
during type I X-ray bursts providing M/R; this is, so far, 
unique to this source. However, these features were not de- 
tected in subsequent observations and thus require further 



confirmation (see Cottam et al. 20081. Recently, Galloway 



et al. ( 2009 1 announced the discovery of nearly periodic os- 



cillations during two type I X-ray bursts in EXO 0748-676 at 
a frequency of 552 Hz. This is at odds with the 45 Hz burst 
oscillations reported by Villarreal & Strohmayer (20041. If 



552 Hz is indeed the neutron star spin frequency the redshift 
measurement is challenged further since rotational Doppler 
broadening tends to lower the contrast unless the orientation 
of the neutron star is favorable, with the spin axis pointing in 
our direction. Using a combination of observations, includ- 
ing t he contested redshift measurement, and theoretical find- 
|2006| argues that the mass in EXO 0748-676 is 



Ozel 



mgs 

2.10±0.28 Mq, which if confirmed would rule out many soft 
equations of state. However, a high spin frequency, such as 
552 Hz, besides challenging the redshift measurement, also 
introduces sizeabl e relativistic corrections to the equations 
and thus result of 1 



Ozel 



(2006). 



Attempts to use narrow He, C and N emission lines that 
are probably originating on the heated inner hemisphere 
of the companion star to get a lower limit on the mass of 



the neutron star are presented in Pearson et al. ( 2006 1 and 



Munoz-Darias et al. (20091. The X-ray eclipse duration for 



EXO 0748—676 in outburst has been measured from X-ray 
studies using the EXOSAT and XMM-Newton satellite (8.3 
minutes, |Parmar et al.|1986| and |Homan et al.|2003| respec- 
tively) . 

Recently, the system returned to quiescence after a 
more than 20 year-long active phase (see Wolff et al. 2008ab, 
Hynes fc Jones|2008| |Torres et al.||2008[ ). Analysing the op- 
tical 7?-band and near-infrared J-band light curves, |Hynes| 
& Jones (2009) find that the residual X-ray luminosity of 
the cooling hot neutron star (see Dcgcnaar et al.|[2009| is 
causing the inner hemisphere of the companion star to be 
heated. 

We here present phase-resolved optical spectroscopy 
and X-ray observations of the source shortly after the end 
of the outburst. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 X-ray observations 

We have observed EXO 0748-676 on Nov. 6, 2008 at 
08:53:10 (UTC) with the European Photon Imaging Cam- 
era (EPIC) on board the XMM-Newton satellite. These data 
were acquired 62 days after the first X-ray observations of 
EXO 0748-676 by |Degenaar et al.|p009| . Hence, the source 
was in quiescence. The observation time was awarded to 



our Directors Discretionary Time proposal. The observation 
identification number is 0560180701. Two of the three imag- 
ing cameras on board XMM-Newton use Metal Oxide Semi- 
conductor (MOS) CCD arrays whereas one uses a pn CCD 
array. Since the MOS cameras receive less than 50 per cent 
of the light due to the presence of the Reflection Grating 
Spectrometers in the light path we present results of the 



pn camera only. The pn camera was employed in Prime- 
FullWindow mode and a medium filter was put in place to 
avoid strong contamination due to bright optical sources in 
the field of view. 

To benefit from the latest calibrations available in De- 
cember 2008 we have reprocessed and analysed the data us- 
ing the SAS software version 8.0.0. Data are excluded if 
the 10-12 keV background count rate for single events in a 
3x3 CCD-pixel region (with pattern 0.0) is higher than 
1.0 counts -1 . The net on-source exposure time is 29.82 ks. 
We have produced bary-centered events lists using the SAS 
tool barycen. 

Next, we have produced a background subtracted source 
light curve with 10 second time bins including events of pho- 
ton energies larger than 0.15 keV. Source photons have been 
extracted from a circular region of 35" radius centered on 
the best known source position. Background photons have 
been extracted from a source-free region of the same size 
and shape on the same CCD. 



2.2 Optical observations 

Long-slit spectra of the optical counterpart to 
EXO 0748-676 were obtained with FORS2, the low 
dispersion spectrograph of ESO's Very Large Telescope. 
Between November 25th and December 26th, 2008, a total 
of 31 spectra were obtained, all with 870 s exposure times 
and using the 1200R grism, covering the wavelength range 
between 5750 A and 7310 A, and a 0'.'7 slit. The detectors 
were read out with 2x2 binning, providing a resolution of 
2.1 A sampled at a dispersion of 0.76 A pix -1 . All spectra 
were obtained during times of good to excellent seeing 
(0'.'48 to 0'.'75). 

The images were corrected for bias and flat-fielded us- 
ing lamp flats. For the sky subtraction we used clean re- 
gions on the slit, fitting it with a polynomial of up to 2nd 
order if that provided the best fit. The spectra were opti- 



mally extracted following the algorithm of Horne|(|l986|) and 



wavelength calibrated using arc lamp exposures taken dur- 
ing daytime. Finally, an approximate flux calibration, not 
correcting for slit losses, was determined from a 22 s expo- 
sure of the spectrophotometric standard LTT3218 taken on 
the first night with a 5" slit. 

We furthermore analysed 15 J?-band acquisition images 
of 30 or 40 s exposure times taken prior to the spectral obser- 
vations. These images were corrected for bias and flat-fielded 
using sky flats. Instrumental magnitudes of stars on each im- 
age were determined through point-spread-function fitting 
using DAOphotll ( |Stetson|p87] . All magnitudes were cal- 
ibrated against the instrumental magnitudes of stars on the 
first image by using some 350 to 400 comparison stars to fit 
for an offset in instrumental magnitude. 

For the photometric calibration of the first image, we 
used aperture photometry of a few bright stars to de- 
termine the aperture correction and determined the pho- 
tometric calibration from 19 photometric standards from 
the PC 2213-006 field ( |Landolt|p92l |Stetson|[2000| . We 
used the standard FORS2 .R-band extinction coefficient of 
0.08 mag per airmass to correct for the difference in airmass 
between the field of EXO 0748-676 and the PC 2213-006 
standard field. Because no color coefficients could be deter- 
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Figure 1. X-ray and optical variability and the radial velocity 
folded on the constant-period orbital ephemeris of |Wolff et al.| 
(2002). Two orbits are shown for clarity, (a:) The X-ray light 
curve showing the eclipse at phase cj> = 0. (b:) R-band magnitudes 
determined from the acquisition images. ( c:) Equivalent width of 
the narrow Ho component, determined from Gaussian fits (see 
§ |3.2| . (d:) The average optical flux from the individual spectra 
between 6000 A and 6250 A. Both in broadband filters and narrow 
lines, the optical flux of the optical counterpart varies with the or- 
bital period, being faintest at the time of the X-ray eclipse. Note 
that variable slit losses are likely present, (e:) Radial velocities 
determined by fitting Gaussians to the Ha profile. Only measure- 
ments between 0.15 < <f> < 0.85 are shown. The grey curve shows 
the best fit circular orbit (v(<f>) = 7 + K cul sin27r</>) to the data. 



mined, we estimate the uncertainty in the photometry as 
the rms in the zeropoint, which is about 0.1 mag. 



3 DATA ANALYSIS 
3.1 X-ray timing 

During the X-ray observations of EXO 0748—676 the source 
shows a constant count rate of 0.55 ± 0.01 counts s except 
for the two distinct eclipses. Fig.[TJi shows the X-ray light- 
curve of the total detected counts (including background). 
The background flux was consistent with being constant at 
0.015 counts s _1 . For the present work we are mainly inter- 
ested in the eclipse duration and the orbital period and mid 
eclipse times. The spectral properties of the source will be 
presented elsewhere. We fit the light curve by a constant 
plus four step functions located at To — |A£ and To + |Ai 
for the first eclipse and T + Pb - §Ai and T + P b + §At 
for the second eclipse. Here To is the mid eclipse time of 
the first eclipse, Pb is the orbital period and At is the 
eclipse duration. The step functions are approximated by 
heaviside functions, defined as f(t) = 5 + 5 tanh(fct), where 
k = —400 hour -1 is used to describe the eclipse ingress 
and k — 400 hour -1 describing the eclipse egress. With 
these values for k, 50 per cent of the flux at the ingress and 
egress times is contained in a single 10 s bin. Fitting for 
To, Pb, At and a constant to describe the flux out of the 
eclipse, we obtain T = 54776. 501663T0. 000068 MJD/TDB, 
P b = 3.8276 ± 0.0018 h and At = 500 ± 7 s. 

We have compared To and Pb with the values predicted 
by the orbital ephemerids of EXO 0748—676 obtained from 
RXTE X-ray timing by | Wolff etH] \20Q2\. Extrapolating 



their ephemerids to orbit n = 54384 gives consistent results 
for the orbital period at the 2a level, while the mid eclipse 
time falls between the predictions of the constant-period 
ephemeris and the quadratic ephemeris using all data. At 
6<T, the prediction from the constant-period ephemeris lies 
closest to our measurement. Because the ephemeris of |Wolff| 
et al. (20021 uses a considerably larger dataset, their results 



are expected to represent the long term timing better, and 
hence, for the remainder of the paper we will use the mid 
eclipse time Tq and orbital period Pb given by the constant- 



(2002) 



period ephemeris of Wolff et al 

We selected in total 800 s worth of data centered on 
the mid-eclipse times for both eclipses in our dataset. Even 
during the X-ray eclipses, EXO 0748—676 is still detectable, 
providing 29 source counts (38 total counts from which 9 are 
expected due to the background) in the 800 s of observation 
that we used. This constitutes approximately 6.6 per cent of 
the out of eclipse 0.3 — lOkeV count rate. 



3.2 Radial velocities and spectral properties 

The individual extracted optical spectra are generally fea- 
tureless except for the region around Ha where all spectra 
have a broad emission component (FWHM of about 26 A), 
while a strong narrow emission line (5 A FWHM) is present 
in spectra taken between orbital phases of approximately 
0.15 < 4> < 0.85. Here 4> — is defined as the center of 
the X-ray eclipse. Fig. [2] shows a trailed spectrum of the Ha 
region, folded on the ephemeris of Wolff et al. ( 2002 I , where 
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the emission clearly follows the orbital signature of a circular 
orbit. The emission appears strongest between the quadra- 
tures and weakest at the time of the X-ray eclipse. 

To measure radial velocities of the Ha components we 
fit the spectral range between 6400 A and 6700 A of each 
spectrum with either a single Gaussian for spectra taken 
between <f) > 0.85 and <j> < 0.15 and two Gaussians for 
spectra taken between 0.15 < <f> < 0.85. Three parameters, 
position, width and height, were fitted for each Gaussian. 
Fitting a circular orbit with v(<j>) = 7 + K eln sin2n<f> to the 
radial velocities of the narrow emission component returns 
= 78.6 ± 3.9 km s" 1 and K em = 333.0 ± 5.0 km s" 1 with 
a reduced xi — 2.4 for 22 degrees of freedom (see Fig.jlJ). 
The errors on the values given above are the formal la uncer- 
tainties and are scaled to give Xt — 1- Because these radial 
velocities are determined from emission lines the measured 
radial velocity semi-amplitude K cul may or may not be equal 
to the true radial velocity semi-amplitude of the companion 
K 2 , see §|4] 

For this fit, the orbital period p, and mid eclipse time 



7 



To were set to those of the ephemeris of Wolff et al. (2002) 



Leaving these parameters free in the fit gives values that 
are fully consistent with their constant-period ephemeris 
extrapolated to n = 54595. The difference in To is only 
1.1 ± 0.7min (0.005 ± 0.003 in phase), while the difference 
in P b is 0.18 ± 0.38 s. 

An average of the 31 individual spectra is shown in 
Fig. [3] where the individual spectra are shifted to zero radial 
velocity using the Ha radial velocity curve. This spectrum 
shows, besides the narrow and broad Ha components, weak 
emission lines of He I A5875, A6678 and A7065. Some tel- 
luric absorption features are present at 6800 A. Zero velocity 
averages of individual spectra with orbital phases between 
0.25 < <f> < 0.75 and between </> > 0.75 and (/> < 0.25 are 
also shown in Fig.[3] Here, the first spectrum predominantly 
shows the hemisphere of the companion facing the neutron 
star, while the latter spectrum shows the opposite hemi- 
sphere. The spectra are significantly different, with the Ha 
and He I emission lines being strongest on the hemisphere 
facing the neutron star. Furthermore, that hemisphere ap- 
pears brighter by about a factor 2. 

Though no absorption features appear to be present 
that may be attributable to the binary companion, we have 
compared the averaged spectrum against synthetic stellar 
spectra by Munari et al. (20051. In the wavelength range 



covered by our observations, the strongest absorption fea- 
tures are the TiO bands around 6180 A and 7050 A for tem- 
peratures below about T e g — 4000 K, while the NalD lines 
are present for temperatures between 4000 K to 7500 K, de- 
creasing in strength with increasing temperature. Between 
4000 K to 5000 K many metallic lines are present in the range 
between 6000 A to 6400 A. Above about 5000 K Ha becomes 
increasingly stronger. 

Unfortunately the region around the NalD lines over- 
laps with the He I A5875 emission line while the Ha emission 
likely swamps the absorption component. As such, both re- 
gions cannot be used. Instead, we use the region between 
6000 A and 6250 A, computing a \ 2 statistic by fitting the 
ratio in flux between the template and the zero velocity av- 
eraged spectrum with a second order polynomial. We find 
that the best fit solution sets a 3ct limit on the temperature 
of Teff > 5000 K, regardless of the surface gravity. This tem- 



perature is the temperature above which the TiO bands, 
disappear. The fit has x? = 1-5 for some 4900 degrees of 
freedom. The absence of absorption lines in the observed 
and synthetic spectra essentially prevent us from obtaining 
a better constraint on the temperature of the companion, 
and measuring its radial velocity curve using photospheric 
absorption lines. 

3.3 Variability 

The spectroscopic as well as the photometric observations 
show that the optical flux varies in both the continuum 
and emission lines. Fig.fTJs shows the P-band magnitudes 
determined from the acquisition images and they confirm 
the light curve published by |Hynes fc J ones (2009); a single 
peaked modulation with the minimum at the time of the X- 
ray eclipse. This trend is also present in both the equivalent 
width of Ha (Fig.JTJ;) and the continuum emission (Fig.[TJi). 
The narrow Ha emission seems to disappear completely near 
the phase of the X-ray eclipse at <fi = 0, while the broadband 
flux varies by a factor of about 2. 



3.4 Doppler tomography 

In order to prepare the spectra for emission line tomogra- 
phy, we first continuum normalised each spectrum using a 
spline-fit and re-binned the data to a constant velocity scale 
sampling of 33.8 kms" 1 pix" 1 . The spectra were then phase- 



binned using the Wolff et al. ( 2002 1 ephemeris. The promi- 



nent orbital variation in the Ha component led us to employ 
the modulation Doppler tomography method described in 



Steeghs ( 2003 \ . This technique includes variable contribu- 



tions modulated on the orbital period and was used to con- 
struct emission line tomograms resolving the Ha and He I 
A6678 emission distributi ons. We used a systemic velocity of 
78 kms" 1 as derived in 



3.2 



although the reconstructions 
are robust against the exact choice of 7. The advantages of 
Doppler tomography are its ability to map complex line pro- 
files using the whole data set at once. Contributions from 
specific sources in the binary can then be identified by their 
known position on the tomogram. 

Reconstructions achieved good fits to the observed line 
profiles, achieving formal \i ~ 1-1 f° r both Ha and He I 
A6678. We show the reconstructions for these lines in Fig. [2] 
The bulk of the emission originated from a region consistent 
with the location of the mass donor star, although a weak 
underlying signature from an accretion disc can also be seen. 
These are responsible for the broad base (disc) and nar- 
row (donor) moving line components previously mentioned. 
Since we are mainly concerned with the contribution from 
the donor, we restrict the plotted velocity intervals to the 
region surrounding the expected location of the donor. 

In Ha, the donor star component is reconstructed as an 
extended feature centred on the Y-axis, covering velocities 
between 200 and 410 kms" 1 . This emission is furthermore 
significantly modulated with an amplitude of 30 per cent, 
peaking near orbital phase <j> — 0.5. Both the distribu- 
tion as well as the modulation of this emission is exactly 
what would be expected for line emission from an irradi- 
ated donor star. However, the emission does not appear to 
be uniformly distributed across the Roche lobe. We mea- 
sured the centroid position of the pronounced peak in our 
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Figure 2. Tomographic reconstructions for Ha (left panel) and He I A6678 (right panel). In each case, the top-left panel shows the 
observed line profiles as a function of binary phase. The top-right panel shows the reconstructed data from the converged maximum- 
entropy solution. Both are displayed on the same gray scale as shown in the wedge. Bottom panels show the corresponding tomograms 
resolving the line emission across the V x — V y plane. Left is the constant contribution to the flux, right the phase-dependent component 
(see[Stccghs 2003]l. The gray scale wedge denotes the fractional amplitude of the variable map in per cent. For reference a Roche lobe 
model and ballistic stream trajectory is plotted for an assumed system mass ratio of q = 0.25 and projected donor star orbital velocity 
of K 2 = 410 kms -1 . 



Doppler tomograms using several methods and found con- 
sistent velocity measurements to within a few kms -1 . For 
example, a Gaussian fit to the Ha map delivers a peak 
of V y — K cln = 413 kms -1 while the highest pixel is at 
410 kms -1 . We also reconstructed Doppler maps using dif- 
ferent values for the systemic velocity to see by how much 
the measured peak position would be affected. Varying the 
systematic velocity 7 by as much as 15 km s -1 shifts the peak 
velocity down by up to 5 kms -1 . The donor star component 
is broadened and weaker for reconstructions with 7 signifi- 
cantly above or below our assumed value of 78 kms -1 . This 
is further confirmation that this is close to the true value for 
the systemic velocity and that its impact on the measured 
K em is small. Given that formal error propagation is not 
possible in maximum-entropy regularised inversions such as 
our Doppler maps, we conservatively assign an uncertainty 
of 5 kms -1 to our measured velocities. This uncertainty en- 
compasses the observed scatter in derived K CYn when com- 
paring different centroid methods as well as the systematic 
contribution due to a possible uncertainty in 7 (as evaluated 
above) . 

For Ha, our measured K cm = 410 ± 5 kms -1 is signifi- 
cantly larger than the K eul values determined in § |3.2| via our 
radial velocity curve fits. For comparison, He I A6678 shows 
a more spot-like distribution, again centred closely to the 
Y-axis and peaking at V y = K em — 345 ± 5 kms -1 . We do 
not have sufficient signal-to-noise in this much weaker line 
to constrain its variability amplitude to a significant degree. 
We thus find that the peak of the Ha emission is found fur- 
ther back on the Roche lobe compared to the He I emission. 
For illustration, we plot the expected shape of the donor 
star Roche lobe in Fig. [2] using a mass ratio of q = 0.25 



and K2 = 410 kms . The various estimates for K BU1 will 
be compared and discussed in the next section. 



4 DISCUSSION AND CONCLUSIONS 

We have obtained X-ray timing and phase-resolved op- 
tical spectra of the eclipsing low-mass X-ray binary 
EXO 0748—676 shortly after the source transitioned into 
quiescence after several decades of outburst activity. Timing 
of the X-ray eclipses show that the orbital period and eclipse 
duration are consistent with the values when the source was 
in outburst. In addition, the orbital phase is consistent with 
extrapolating the constant period ephemeris of | Wolff et al.| 
p002| . 

The X-ray eclipses in EXO 0748—676 are near total, 
however, about 6 to 7 per cent of light remains during the 



eclipses. Bonnet-Bidaud et al. (2001) find that in outburst 



part of this residual flux can be attributed to fluorescent 



thermal halo emission. Homan et al. ( 2003 1 find that in out- 



burst about 4 per cent of the 2 — 6 keV continuum emission 
is detected during eclipse. They attribute this to light scat- 
tered into our line of sight by a corona. Given that they 
find that the fraction of light during the eclipse is larger for 
lower energies, our value seems consistent with that found 
in outburst. This implies that the structure responsible for 
the scattering is still present, or still replenished even when 
the source is in quiescence. 

Using our medium resolution VLT/FORS2 spectra of 
EXO 0748-676 we find emission lines of Ha and He I. They 
display the orbital motion expected if the lines originate on 
the binary companion. The optical spectra show consider- 



6 Bassa et al. 



7 

eg 
I 

E 

o 

o> ° 

L. 

CD 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - 

He 1 


1 1 1 1 1 1 1 1 1 1 1 1 

He 1 He 1 


9 DIB 1 H « ® 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - 

O.25<0<O.75 
(6<0.25, 0>O.75 


1 1 1 1 1 1 1 1 1 1 1 1 


■ i i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i 



6000 



6250 



6500 
Wavelength (A) 



6750 



7000 



Figure 3. (top:) The spectrum of the counterpart to EXO 0748— 676. This spectrum is the average of the 31 individual spectra, each 
shifted to zero velocity using the orbital fit to the narrow Ha emission. The broad and narrow Ho components, as well as three He I 
emission lines are visible. Regions of telluric absorption and diffuse interstellar bands are indicated, (bottom:) Zero velocity averages 
of the individual spectra with orbital phases between 0.25 < <j> < 0.75, basically showing the spectrum from the hemisphere of the 
companion that is facing the neutron star, and between <j> > 0.75 and <f> < 0.25, showing the spectrum of the opposite hemisphere. It 
shows that the emission lines are strongest on the hemisphere facing the neutron star, and that that hemisphere is brighter by about a 
factor 2. 



able flux variations in both the continuum and the emission 
lines. The variations display a single peak over the orbital 
period and we attribute this to irradiation of the hemisphere 
of the companion facing the neutron star. |Hynes fe Jones| 
(2009) also presented evidence for an irradiated donor star 



on the basis of their optical and near-infrared photometric 
observations. 

This irradiation causes three effects that will skew the 
measurement of the radial velocity semi-amplitude. First, 
due to the irradiation the centre of light is offset from the 
centre of mass of the companion. The former being on the 
hemisphere facing the neutron star, the measured radial ve- 
locity semi-amplitudes are smaller than the actual orbital 
velocity of the companion. To correct for this difference one 
can apply a so called 'Tf-correction' ( |Wade fc Home 1988[ ). 
This 'iC-correction' can be estimated assuming values for the 
mass-ratio q and the opening angle of the accretion disk. The 
latter is needed to correct for the partial shielding of the bi- 
nary companion by the accretion disc from the irradiating 
flux (cf. Munoz-Darias et al.|2005 l. Second, the irradiation 
may induce gas motion on the surface of the mass donor 
star, resulting in apparent shifts in orbital phase of the ra- 
dial velocity curve (Kirbiyik 1982 cf. Jonker et al. 20051. 
Third, due to the difference in flux observed from the ir- 
radiated and non-irradiated hemispheres of the companion 
the rotational broadening profiles are not symmetric. As a 
result observed line profiles can become dependent on or- 
bital phase and inclination ( Jonker et al.|2003 Jonker et al 



|2005| ) . This effect is obvious when considering the extreme 
case where all the flux originates from the irradiated hemi- 
sphere. Assuming synchronous rotation, at phase cj> = 0.25 
the companion is moving away from the observer at its max- 
imal observed orbital velocity, however, the observed irradi- 
ated hemisphere is rotating towards the observer lowering 
the observed semi-amplitude of the radial velocity curve. 

How much the measurement of the radial velocity semi- 
amplitude is influenced depends therefore on whether the 
semi-amplitude is determined by fitting Gaussians to the 
emission line profiles or by fitting a two dimensional Gaus- 
sian to the spot in a Doppler tomogram. We modelled the 
changes in the line velocities by fitting a double Gaussian, 
one narrow, one significantly broader, to the Ha emission 
profiles. We find that the velocities of the narrow component 
can be well described by fitting a sinusoid with a constant. 
The semi-amplitude is 338.8 ± 4.8 km s -1 and the constant 
representing the systemic velocity 7 = 78.5 ± 3.4kms _1 
when restricting the fit to the velocities over the orbital 
phases 0.25 < cj> < 0.75 when the narrow component of 
the emission line is strongest. 

Using Doppler tomography of the Ha emission line with 
7 = 78kms _1 as input we find that the centre of the Ha 
emission spot in the Doppler map falls at K em = 410 ± 
5kms _1 . This value for the radial velocity semi-amplitude 
measured from Doppler tomography is significantly larger 
than that obtained fitting a double Gaussian to the emis- 
sion peak. The reason for this difference lies in the fact that 
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the Doppler tomography method does not require symmet- 
ric line profiles whereas fitting Gaussians obviously does. 
The Doppler maps do show that Ha emission is complex 
in shape and thus that the single Gaussian model for the 
narrow component is likely too simplistic to reproduce the 
observed distribution. The third effect mentioned above ren- 
ders the line profiles asymmetric, affecting the determined 
velocities for the 'double Gaussian' method predominantly. 
The Doppler tomography will, however, suffer less from this 
effect yielding velocities closer to the true centre-of-mass ve- 
locities. Note, however, that also the Doppler tomography 
velocities follow the centre of light and hence they also re- 
quire an unknown './^-correction'. 

For the only other emission line strong enough for a 
Doppler map study, the He I A6678 line, we obtain K cul — 
345 ± 5kms _1 . There are a number of effects controlling 
the distribution of line emission across the Roche lobe, and 
hence defining the location of the centre of light. Parts of 
the lobe may be shielded from irradiation due to a vertically 
extended disc, and the optical thickness of this will be de- 
pendent on wavelength. The incident angle of the incoming 
radiation is also much larger near the terminator compared 
to the region near the inner Lagrangion point (LI). Different 
emission lines will thus show a different distribution reflect- 
ing the varying local conditions of the gas, the ionisation 
potential of the element involved as well as the geometry of 
the incoming radiation field. This typically leads to higher- 
ionisation potential lines to be formed closer to the LI point 
(e.g. |Harlaftis||l999l |Unda-Sanzana et al. 12006] ) . 

Our finding that He I emission peaks closer to the LI 
point compared to Ha thus agrees with this picture. It shows 
that the Balmer lines can still be excited near the termina- 
tor, but at the same time also suffer from more effective 
shielding near the front of the lobe because of the increased 
optical depth of the accretion flow due to neutral Hydrogen. 

The measured radial velocity semi-amplitudes are 
higher and carry smaller statistical errors than those found 
during outburst by Munoz-Darias et al. (20091. Besides a 
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higher radial velocity semi-amplitude, we also find a higher 
systemic velocity, although the values are consistent within 



the 3a errors. Munoz-Darias et al. (20091 observed the sys- 



tem in outburst using emission lines from He II whereas our 
observations are done in quiescence and using Ha. The dif- 
ferences in X-ray flux between outburst and quiescence and 
the different element, may explain the difference in mea- 
sured radial velocity semi-amplitude in our data compared 
to 



Munoz-Darias et al. (20091 as explained above 



Line emission from the secondary can be expected to lie 
anywhere between its LI point and the terminator. Rather 
than make specific assumptions about the irradiation geom- 
etry, we will instead use the peak of the Ha emission from 
the Doppler map as a determination of the radial velocity 
semi-amplitude, since it provides the strictest lower limit on 
the true radial velocity semi-amplitude of the mass donor, 
i.e. K2 > K cul . Hence, with K cm = 410 ± 5kms~ 1 , we find 
that K 2 > 405 km s" 1 . 

So far, we have used the fact that the true K2 is larger 
than the observed K cul , where the 'Tf-correction' depends on 
the distribution of the emission across the Roche lobe. Thus 
the lines with the largest K aul must be formed closest to the 
terminator, whereas lines with lower K eul are more towards 
LI. An additional limit can thus be set by demanding that 
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Figure 4. Constraints on the neutron star mass in 
EXO 0748—676. Under the assumption that the companion is fill- 
ing its Roche lobe, the eclipse measurement of the eclipse duration 
gives a relation between the mass ratio q = M2/M1 and the in- 
clination of the orbit (Home 1985). Furthermore, the mass ratio 
relates to the rotational velocity v sin i and the radial velocity 
semi-amplitude The dashed lines show curves for specific ro- 
tational velocities. Finally, the neutron star mass function relates 
the neutron star mass to the mass ratio and the radial veloc- 
ity semi-amplitude, giving the solid lines for specific neutron star 
masses. The horizontal grey area is excluded by our lower limit 
on the radial velocity semi-amplitude of K cul = 410 ± Skms^ 1 , 
setting K2 > 405 km s -1 , while the diagonal grey area is ex- 
cluded in order to satisfy that all emission originates between 
the radial velocity of LI and the center-of-mass of the companion 
(Vli < A' cm < X2). These constraints dictate that the neutron 
star is heavier than Mi > 1.27Mq. The hashed area is the mass 
range of Mi = 2.10 ± 0.28 M favoured by the work of 
(2006) ■ 



Ozel 



Vli < -K"em. |Mufioz-Darias et al.| ( [2"00 9i) reported K cul values 
that are somewhat lower than ours and we used their lower 
limit on K arn = 300kms _1 to calculate which mass ratios 
are ruled out because the velocity of the LI point would be 
larger than 300 km s" 1 , which would be inconsistent with the 
donor star components observed during outburst. This leads 
to a lower limit on the system mass ratio q for a given K2 
which is indicated by the shaded region in the upper right 
of Fig. [4] From this we see that the range of measured donor 
star velocities that all must satisfy Vli < K cln < K2 rule out 
neutron star masses below 1.27M@, so that Mi > 1.27Mq. 
For the canonical neutron star mass of Mi = 1.4 Mq our 
observations constrain the mass ratio to a narrow range of 
0.075 < q < 0.105 and the companion mass to 0.11 Mq < 
M 2 < 0.15 M Q . 

Narrow Ha emission originating from the secondary 
star has also been observed in other quiescent low mass X- 
ray binaries. When measured after subtracting the underly- 
ing Ha absorption line from the secondary and accounting 
the veiling from the accretion disc, the equivalent width in 
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these systems is of 2-3 A and comparable to that observed 
in rapidly rotating chromospherically active stars (see e.g. 
Casares et al.|1997||Torres et al.|2002||2004| ). In contrast, the 
narrow Ha component in EXO 0748—676 reaches an equiv- 
alent width of tenths of A (Fig.JT]). These high values of the 
equivalent width can be explained as due mostly to X-ray 
irradiation of the secondary. Stellar rotation in G-M dwarfs 
cannot produce such a strength because of saturation effects 



of stellar activity with rotation (see e.g. |Soderblom et al 
1993{|Stauffer et al.||1997 l. 

Fitting theoretical stellar template spectra to the ob- 
served spectrum averaged in the frame of the compan- 
ion star, we constrain the temperature of the inner face 
of the mass donor star to be larger than 5000 K. Stars 
of such a temperature do not show significant absorption 
lines in the spectral range covered by our observations. 
Such temperatures are unexpected for the orbital period 
of EXO 0748—676. A Roche lobe filling main sequence star 
would have a mass corresponding to that of an M2 star and 
thus a temperature of approximately 3500 K, much lower 
than the lower limit derived above. Besides evidence from 
the equivalent widths of the above mentioned emission lines 
and the comparison of theoretical stellar templates, we find 
that the continuum varies in accord with a heated inner 
hemisphere of the mass donor star (see Fig.JTJl), hence, the 
variability in magnitude over the orbit cannot be explained 
fully by the change in equivalent width of the emission lines. 

Repeat observations targeting a different wavelength 
range showing deep stellar absorption lines from the mass 
donor star even in the presence of heating due to the hot neu- 
tron star have been requested. Furthermore, the hot neutron 



star probably causing the irradiating flux (Hynes & Jones 



2009) could cool significantly on timescales of a year reduc- 
ing the irradiating flux (Dcgenaar et al.|2009[ ). 
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